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ABSTRACT 



We have imaged 45 quasars from the Sloan Digital Sky Survey (SDSS) with 
redshifts 1.85 < z < 4.26 in JHK S with the KPNO SQIID imager. By com- 
bining these data with optical magnitudes from the SDSS we have computed 
the restframe optical spectral indices of this sample and investigate their rela- 
tion to quasar redshift. We find a mean spectral index of (a ) = —0.55 ± 0.42 
with a large spread in values. We also find possible evolution of the form a Q = 
(0.148 ± 0.068)^ - (0.964 ± 0.200) in the luminosity range -28.0 < Mj < -26.5. 
Such evolution suggests changes in the accretion process in quasars with time 
and is shown to have an effect on computed quasar luminosity functions. 

Subject headings: Quasars and Active Galactic Nuclei: evolution 
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Introduction 



Th e recent publicat i on o: 



(SPSS; 



Schneider et al. 



Croom et al. 



large samples of quasars from the Sloan Digital Sky Survey 



2003) and the Two Degree Field QSO Redshift Survey (2QZ; 



2004 ) allows the study of the evolution of the quasar population spanning 



long lookback times with large, homogenous c atalogs of object s and t he characterization 



Richards et al. 



(j2006al ) QLF formed from 



of the quasar luminosity function (QLF). The 
the SDSS Third Data Release Quasar Catalog (DR3QLF) is defined from 15,343 quasars 
over 1622 deg 2 , spanning redshifts from z = to 5 and reaching luminosities down to 
Mi < —23 in their lowest redshift bins. They find a peak in type 1 quasar activity between 
z = 2.2 and 2.8 and a general flattening of the slope of the bright-end QLF with increasing 



redshift. iHopkins et al. I (j2007f l compute a bolmetric QLF, combining many quasars surveys, 
including the DR3QLF and the 2QZ. They find a peak in the quasar luminosity density at 
z = 2.15 and an evolving QLF bright-end slope that becomes flatter at redshifts z > 3. 

One major difficulty in characterizing the evolution in quasar space densities is finding 
an appropriate way to calculate the absolute magnitudes of a survey sample and compare 
these magnitudes with other surveys, or even to objects within the same survey if it spans 
a large redshift range, i.e., how to calculate the appropriate K cor rection. Traditionally, 



absolute magnitud es were corrected back to the restframe S-band (ISchmidt fc Green 



Boyle et al. 



1983 



19871 ). as most early surveys for quasars included the measurement of flux 
around 4400A, by assuming a power-law form for quasar spectra redward of Lyman-a of 
the form j v oc v a and assuming an average spectral index, a, for the survey sample. As 
surveys move to higher redshifts, beyond z = 3, the observed S-band ceases to sample the 
quasar continuum effectively, as the Lyman-a line moves into and redward of that filter. 
In the recent past, research groups have adopted several methods for computing absolute 
magnitudes and have referenced their magnitudes to different restframe wavelengths. For 
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instance. ISchmidt. Schneider fc Gunn I (119951 ) and lKennefick et al. I (|19951) adopt an average 



spectr al index of a = —0.5 in their computations of Mb, while 



Warren. Hewett fc Osmer 



( 119941 ) avoid the adoption of a spectral index by computing the QLF as a function of 



M, 



C(1216)? 



the flux in the continuum under the Lyman-a line, fr om their direct measu rements 



Croom et al. 



of the flux at that point in their survey spectra. More recently, 
their QLF from the 2QZ at 0.4 < z < 2.1 in terms of M&j, and iRichards et al. 



( 200411 prese nt 



(12006af l 



assume an average spectral index, but correct z-band magnitudes to z — 2. The correction 
to z = 2 is prompted by the desire to minimize the effe cts of extrapolati ng the assumed 
quasar powerlaw over large wavelengths, as suggested by 



Wisotzki 



(120001 ) . 



Such differences in the methods for computing absolute mag nitudes has led to some 



1993). 



discrepancies in computed quasar space densities in past surveys (iKennefick et al. 
The initial goals of this near-infrared (NIR) imaging program were to measure the restframe 
.B-band flux for a set of quasars in several redshift ranges, to compare the differences in 
the absolute magnitudes, Mb, as computed from extrapolations of flux values at lower 
wavelengths to flux values measured directly, and to see how this changed with redshift. 
This is equivalent to measuring the restframe optical spectral index of the quasars from 
optical and NIR photometry. In this paper we report the initial results of a program 
to measure the restframe optical spectral index, a Q , of a subset of SDSS quasars at 
1.85 < z < 4.26 using their reported optical magnitudes and NIR photometry obtained at 
the KPNO 2.1m telescope using the Simultaneous Quad Infrared Imaging Device (SQIID). 
In §2 we present the program design. The data used in the project, both archived and new 
observations, are described in §3. In §4 we give the program results, including computations 
of spectral indices and the implications for the resulting absolute magnitudes. We discuss 
our results further in 55. 
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2. Program Design 

The initial goal of this observing program was to measure the restframe 5-band flux 
for a sample of quasars at a range of redshifts in order to essentially bypass the need to 
extrapolate a power-law over large wavelength ranges when applying a K correction in the 
computation of quasar absolute magnitudes, Mb- The central wavelength of the J-band is 
1.27 microns, which corresponds to the restframe 5-band (central wavelength of 4400 A) at 
z = 1.88. Likewise, if-band corresponds to the restframe 5-band at z — 2.80, and i^-band 
at z — 4.06. Figure 1 shows the NIR JHK S filter curves, along with the "redshifted" 
Johnson 5-band at these redshifts, that is, the portion of the spectrum you would like to 
measure to sample the restframe 5-band at the given redshift. For this project, we targeted 
45 quasars chosen from the SDSS quasar catalogs for NIR imaging centered around the 
above redshift ranges. We combine this NIR photometry with reported optical magnitudes 
to compute restframe optical spectral indices for the sample. Programatically, this can 
be achieved by fitting the photometric data to measure the slope of the spectral energy 
distribution, a, as dicussed further below. 



Richards et al. 



(l2006al ) has reported the DR3QLF computed from SDSS quasars as 
a function of Mj corrected to a redshift of z = 2, and the selection of the S-band for 
computations of the QLF has become less common as quasar surveys move to higher 
redshift. Here, we compute the optical spectral index, a , for a subset of the SDSS quasars 
using optical and NIR photometry. We then explore possible correlations of a Q with redshift 
or luminosity and how that might affect the characterization of the QLF and its evolution. 
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3. Observations and Data Reductions 
3.1. SQIID Near-Infrared Imaging 

Forty-five quasar s from the SDSS Third Data Release Quasar Catalog (DR3Q; 



Schneider et al. 



20051 ) were imaged with SQIID on the KPNO 2.1m telescope during 2005 
March 22-23. The SQIID Infrared Camera uses individual 512 x 512 pixel quadrants of 
ALADDIN InSb arrays, with a pixel scale of 0.69"pixel _1 at the KPNO 2.1m. The effective 
field of view is 304" x 317". We configured SQIID to take images in the JHK S bands, which 
it acquires simultaneously by splitting the beam with a series of dichroics before sending 
the beam to separate NIR cameras optimized to perform over a limited wavelength range. 

The observing strategy involved imaging each target quasar for a total of 600s, with 
five pointings of 120s, each the sum of fifteen 8s coadds, offset by roughly 45" between 
pointings to improve background subtraction. Conditions were nonphotometric, with light 
cirrus over the course of the run. Seeing ranged from 1.3" to 1.9" over the two nights of the 
observing run. 

The images were processed using the IRAfQ UPSQIirJ] package. Dark frames were 
acquired for each combination of coadd and exposure time utilized during the run. A global 
flatfield for each night and each filter was constructed from the object frames using the 
USQFLAT routine, configured to subtract the dark current and then to perform a median 
combine of the frames. Image processing was performed using the MOVPROC routine to 
generate and subtract a moving sky image created from the 6 frames obtained closest to the 
object frame in time and to correct for pixel to pixel sensitivity differences by dividing by 



X IRAF is distributed by the National Optical Astronomy Observatory, which is operated 
by the AURA, Inc under cooperative agreement with the NSF. 
2 http: / /www.noao.edu/kpno/sqiid/ 
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the appropriate flatfield. The final image of each quasar in each band was constructed by 
registering the five pointings using the USQMOS routines and combining the frames using 
the NIRCOMBINE routine. 

Object detection and measurement was carried out using the Source ExtractoiP 
software. The observations were made under fair but nonphotometric conditions. 
Calibration was p erformed through the use of stellar Two Micron All Sky Survey (2MASS; 



Cutri et al. 



20031 ) sources in the object frames. The 2MASS images and catalogs were 
accessed through the National Virtual Observatorjo (NVO) Op en SkyQuery and Da taScope 



20001). The 



Query tools and manipulated using the Aladin multiview tool ( IBonnarel et al. 
JHK S magnitudes reported in Table 1 and used in the following analysis were derived 
from the flux inside circular apertures of radius 3.5". The photometric zeropoints for 
the frames were calculated using stellar 2MASS sources in the frame of the quasar with 
measured JHKs aperture magnitudes. Typically five to seven calibration stars were used 
with magnitudes ranging from 14 to 16 in J and H, and from 14 to 15 in K s . The 2MASS 
aperture magnitudes were measured in a 4" radius and "curve-of-growth" corrected out to 
an "infinite" aperturqj. Therefore, the magnitudes given in Table 1 are effectively aperture 
corrected by using these corrected 2MASS sources as standards. Near-infrared magnitudes, 
uncorrected for Galactic extinction, for the quasars are given in Table 1 and Figure 2, 
along with their associated errors, which include both measured photometric errors for the 
quasars and the errors introduced by the calibration process. The quasars in the subsample 
observed with SQIID are 1 to 2 magnitudes fainter than the SDSS quasars detected by 
2MASS. The relatively few numbers of quasars in the highest redshift bin (~ 4.06) meant 



3 http:/ /terapix. iap.fr/ 
4 http: / /www.us- vo.org/ 

5 http: / / www.ipac.caltech.edu/2mass/releases/allsky/doc/ 
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that we had to choose fainter objects over a slightly wider redshift range than in the two 
lower redshift bins. 



3.2. The SDSS Third Data Release Quasar Catalog 



The SDSS uses a multicolor technique to select quasar candidates in the optical ugriz 
bands over the redshift range 0.08 < z < 5.41. The DR3Q contains 46,420 quasars with 
luminosities Mj > —2 2. Initially, we chose 6 quasars from the SDSS Early Data Release 



quasar catalog (EDR; 



Schneider et al. 



20021 ) for NIR imaging. Quasars were selected from 
the EDR catalog randomly from those objects with r.a. between 9 h and 16 h and as close 
to the target redshifts as possible. Since the EDR contains more objects at lower redshifts, 
there is less spread in the sample redshifts at z — 1.88, increasing slightly at z — 2.80 and 
more pronounced at z = 4.06 (see Figures 2 and 3). The mean redshifts for the three bins 
are z = 1.88, 2.82, and 4.03. We obtained NIR imaging for 45 of these objects, chosen 
to sample our redshift ranges and with RA's close to the local sidereal time during the 
observations. All of the imaged quasars are contained in the SDSS DR3Q catalog, and the 
optical photometry reported in Table 1 and used in the data analysis were taken from this 
later catalog. 

The SPSS DR3Q has been matched to the 2MASS All-Sky Data Release Point Source 



Catalog (ICutri et al. 



20031 ) and the 2MASS JHK S magnitudes and errors are reported for 
those SDSS quasars with 2MASS detections (columns 25-30 of the DR3Q.) We have used 
the 6192 quasars with 2MASS detections as a bright comparison sample to our fainter (l m 
to 1.5 m ) SQIID sample (Figures 2 and 3). 
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Results 



4.1. Colors 



The magnitudes were corrected for Galactic extincti on using the value of A u for each 



object reported in the SDSS DR3Q catalog, column 15 ( 



values were taken from the maps of 



Schlegel et al. 



Schneider et al. 



20051 ). These 



(119981 ) which assumes an extinction 



to reddening ratio Ry = 3.1. Extinctions in the Sloan griz bands are then 0.736, 0.534, 
0.405, and 0.287 times A n , respe c tively . The value of A/E(B — V) for the Sloan w-band is 



5.155; consulting 



Schlegel et al. 



(119981 ) for the UKIRT JHK values gives 0.902, 0.576, and 
0.367, respectively. This gives values for the extinction in these bands of Aj = 0A75A U , 
An = 0.112Au, and A K = 0.071v4 u . Colors computed from these extinction corrected 
magnitudes for the SQIID sample and the SDSS quasars detected by 2MASS are shown 
in Figure 3. The SQIID sample spans the same range in quasar colors as do the 2MASS 
detected quasars. 

Also shown in Figure 3 are expected colors for quasars computed using synthetic 
quasar spectra generated with power law spectral indices and averaged over several different 
realizations of the Lyman-a forest. The colors were computed by passing the quasar spectra 
through SDSS and SQIID filters and calibrated with a comparison Vega spectrum from the 



Space Telescope Science Institute CALSPECc 
then converted to the SDSS AB system using 



site. The uqriz Vega based magnitudes were 



Holberg fc Bergeron 



(|2006|). The dotted line 



in Figure 3 represents quasars with a = —0.5. While the quasar colors do cluster around 
this line, the colors span a conside rable range, large r than the errors in their colors. This 



is consistent with results found by 



Pentericci et al. 



(120031 ) . who found a similar spread in 



optical/NIR colors in a sample of 45 SDSS quasars at 3.6 < z < 5.03. 



3 http://www.stsci.edu/hst/cdbs/calspec.html 
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4.2. Optical Spectral Indices 



In order to compute the restframe optical continuum slopes of the quasar spectral 
energy distributions, we have converted our NIR, Vega based photometry to the AB system. 
To accomplish this, we first construc ted a NIR spectru m for Vega based on the absolute 
flux calibrations for Vega reported by 



Megessier 



(119951 ) (see their Table 4,) constraining the 



blue end to have f v = 2446JV at A = 9000A, as given in their formula for the blackbody 
fit to Vega (see their page 776.) We then compared the flu x from this spectrum with a flat 



source with /„ = 3631Jy, the zeropoint of the AB system (jOke fc Gunn 



1983|), for each 



2MASS and SQIID JHK S filter. The offsets for the 2MASS filter set are: J AB = J + 0.84, 
H AB = H + 1.33, and K AB = K s + 1.79. For the SQIID filter set they are: J AB = J + 0.88, 
H AB = H + 1.35, and K AB = K s + 1.75. 

If we assume that the optical fl ux for a quasar is giv en by f v = u a ° , then we can rewrite 



the formula for the AB magnitude (jOke fc Gunn 



19831 ) 



AB = -2.5 log 10 f u - 48.6 



in the following form: 

AB = -2.5a o log 10 z/-48.6, (2) 
allowing us to compute a Q as the slope of a fit to straight line in AB magnitude vs. log y. 



19921 ) 



This was accomplished by using the Numerical Recipes routine FIT (IPress et al. 
Spectral indices for each quasar observed with SQIID are given in Table 2. The mean 
spectral index for the sample is (a a ) = —0.55 ± 0.42. However, there appears to be a change 
in the mean slope with redshift in the SQIID sample, with quasars at lower redshifts having 
a steeper slope. In Figure 4, AB magnitudes transformed to the quasar restframe and 
normalized to 20.0 in the bluest band completely redward of Lyman-a are given along with 
a line showing the mean slope for that redshift range. The mean spectral index for each of 
the redshift bins is: (a )i. 88 = -0.71 ±0.43, (« ) 2 . 82 = -0.49 ±0.40, (a o ) 4 . 3 = -0.40 ±0.39 
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(see also Table 3). These spectral indices were calculated using all available passbands. 
For consistency, if we only use those passbands available in all three redshift bins (from 
restframe ~ 1500A out to ~ 4500A), the computed spectral indices change very little, with 
means of: (a )i. 88 = -0.68, (a ) 2 .82 = -0.50, (a o ) 4 .03 = -0.40. 



4.3. K corrections and Absolute Magnitudes 



In order to compute a luminosity function for quasars, the absolute magnitude of each 
discovered object must be calculated. This can be a bolometric luminosity, the luminosity 
at a point in the continuum, or the flux in a passband. However, since the spectra of 
quasars are redshifted and even an individual survey can span a large range in redshift, 
the calculation must include a conversion from an observed band to an emitted band, 
his conversion from observed flux to restframe flux is referred to as the K correction 



(IHumason. Mayall. fc Sandage 



19561 ). defined as the "technical effect that occurs when a 



continuous energy distribution F (X) is redshifted through fixed spectral-response bands of a 



detector" (IQke fc Sandage 



19681 ). If the spectral energy distribution (SED) is not flat, this 



will include both the effect of detecting light from a region of the emitted spectrum shifted 
from the effective wavelength of the detector and the effective squeezing of the detector 
bandpass in the emitted frame. 



Following the formalism of iHogg et al. I (120021 ). the K correction is defined as 



M em itted — m observed ~ DM — K(z) 



where DM is the distance modulus, given as 

DM = 5 log 10 



Dl 
lOpc 



(3) 



(4) 



K(z) is the K correction, which depends on the SED of the o bserved objec t (see 



2002 



eq. 8), and Dl is the luminosity distance. Following 



Hogg 



Hogg et al. 



(120001 ) and assuming 
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Om + ^a = 1, the luminosity distance is given by 



c(l + z) 

d L (z]tt M MA,Ho) = — x 

-no 



[(1 + z') 2 (l + n M z') - z'(2 + z'W^dz'. (5) 







We a ssume H =70km s Mpc , Vt M = 0.3, and £l\ = 0.7 throughout (jSpergel et al. 



20071 ) . 

For a power-law SED of the form f v oc v a we can write 

K(z) = -2.5a log(l + z) - 2.5 log(l + z) (6) 

where the first term transforms from the emitted band to the observed band at z=0 and the 
second term corrects for the effective narrowing of the passband when observing redshifted 
objects. Quasar spectra also contain various broad emission lines, and their contribution to 
the K correction will be addressed below. 

Table 2 (column 4) lists the absolute magnitude, Mj( 2= o), for each quasar as reported in 
the SDSS DR3Q catalog. These values were computed by correcting the SDSS i magnitudes 
for Galactic extinction, assuming H =70km s _1 Mpc" 1 , Qm — 0.3, and Q\ = 0.7, and 
applying a standard continuum K correction by assuming a spectral index of a = —0.5 and 
correcting to z = 0. 

The DR3QLF is computed as a function of Mj, but K corrected to z = 2, near the 
peak in the quasar distribution. They also correct for emission lines by convolving a 
composite spectrum created from 16,713 SDSS quasars with the continuum subtracted with 
the SDSS filters. Their values are reported in Table 2 (column 5) as M,^ z= 2). They compute 
a combined K correction, assuming a spe ctral index of a = —0.5 (see their Table 4), which 



we will refer to here as K2. We adopt the lRichards et al. 



(j2006al ) methodology but employ 



the spectral index computed for our SQIID subsample in computing the Mj. This gives an 
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expression for K(z): 

K{z) = K 2 - 2.5(a - a fix ) log 10 (l + z) (7) 
where af ix = —0.5. Therefore M^ ao ^ is given as 

-K 2 + 2.5(a o + 0.5)log 10 (l + z). (8) 

The values of a Q and Mj( Qo ) for the SQIID sample are listed in Table 2 (columns 6 and 
8), and the M^ ao ) are plotted vs. redshift in Figure 5. Because our imaged z ~ 4 quasars are 
fainter than those in the two lower redshift ranges, the program samples quasars of similar 
luminosity at all three redshifts, unlike the 2MASS crossmatches, which have a distinct 
dependence of luminosity on redshift, as would be expected for a flux limited survey. 



M 



i(cto) 



Tflj 



5 log 



10 



Dl 
lOpc 



4.4. Correlations 



In order to investigate possible correlations of a Q with redshift or 



used the Astronomy SURVival Analysis tools (ASURV; 



La Valley et al. 



uminosity, we 



19921 ). accessed 



through the IRAF STSD AS package. Linear re gressions were performed using both the EM 



(estimate and maximize; 



Dempster et al. 1119771 ) and the Buckley- James (IBuckley fc James 
19791 ) algorithms resulting in the following expressions for a in terms of z and M^-v. 



a Q (z) = (0.148 ± 0.068)^- (0.964 ± 0.200) 



(9) 



a (M i(ao) ) = (0.168 ± 0.083)M i(Qo) + (4.084 ± 2.305) (10) 

a (z,M i{ao) ) = (0.249 ± 0.065)^ + (0.295 ± 0.080)M i(Qo) + (6.891 ± 2.137). (11) 

The expressions for a (z) and a (Mj( Q , o )) are plotted along with the data in Figure 6. 
We have also computed the gener alized Kendall's tau co rrelation coefficent between these 



variables using the BHK method (jLaValley et al. 



19921 ). For the a a and z relation, we 
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find r = 0.40 with a 95% probability of a correlation. For a a and M^q, \, we find r = 0.30 
with an 85% probability of a correlation. We note, however, that the SQIID sample was 
not chosen to sample a broad range of absolute magnitude. Instead fainter sources were 
necessarily chosen at higher redshifts because fewer bright quasars are present at these very 
high redshifts in the SDSS sample due to the relatively fewer numbers of quasars at these 
epochs. We have also included an expression for a (z, Mj( ao )) and show the projections 
of the residuals of this fit along with the residuals from the expressions for a Q (z) and 
a (Mj( ao )) in Figure 7. The residuals in the two-dimensional expressions (top panels) are 
similar to those for the three-dimensional expression (lower panels). 

The distributions of the spectral indices for the SQIID and 2MASS detected samples 
are shown in Figure 8 (left panel). Also shown are the separate distributions of the SQIID 
subsamples at z ~ 1.88, 2.82, and 4.03 (right panel). The ot Q of the 2MASS detected 
SDSS quasars were computed using the same method as for the SQIID sample, using 
the optical bands completely redward of Lyman-a and the available 2MASS data from 
the DR3Q. Sample statistics are given in Table 3. The median of the complete SQIID 
sample distribution is flatter than the mean with a value of —0.47. The mean is shifted 
to steeper values by the presence of a red tail to the distribution. The median of the 
2MASS sample distribution is —0.57. The mean is again shifted redward by a red tail to 
(a Q ) = —0.70 ± 0.53. The distributions of th e two samples look very similar. However, 



performance of a Kolmogorov-Smirnov (K-S; 



Babu fc Feigelson 



1996 



Press et al. 



1992) 



test on the two data sets shows only a 38% probability that the two distributions are drawn 



from the same parent distribution. The differences could be d ue to dependence o 



a on 



luminosity, as has been suggested by this work and others (e.g. ISteffen et al. II2006I ). as the 
SQIID and 2MASS samples cover different areas of luminosity/redshift space (Figure 5.) If 
we consider just those quasars from the 2MASS sample with —28 < Mi( ao \ < —26.5 (where 
the bulk of the SQIID sample lies), then the mean of the distribution is (a Q ) = —0.54 ±0.37 
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and the median is -0.47, almost identical to the SQIID sample, even though the bulk of 
these quasars are at redshifts between 0.5 and 2 (Table 3). 

As is clear from Figure 5, the 2MASS detected sample has a strong correlation of 
luminosity with redshift. While this is also true, to a lesser degree, for the SQIID sample, 
in order to select a statistically significant sample at high redshift and stay close to our 
desired redshift ranges, we had to observe fainter candidates at the higher redshift bins, 
as can be seen in Figure 2. This resulted in the SQIID subsamples having very similar 
mean luminosities of M i{cto) = -27.15, -27.89, and -27.87 at z ~ 1.88, 2.82, and 4.03, 
respectively. It is therefore not surprising that we see stronger evidence in favor of a 
correlation of a Q with z than with luminosity, as our survey was designed to explore the 
relation with redshift. We do note, however, that if we limit the 2MASS detected sample 
to have similar absolute magnitudes to our sample, then the statistics for a a become very 
similar (Table 3), even though the mean reshifts are 1.44 and 2.81. What we conclude is 
that there is evidence for evolution of a Q with both z and luminosity, but that more data at 
a broader range of luminosities is needed to fully characterize its form. 

Computing the difference i n the absolute magnitud es in Table 2, one computed with 
an average spectral index as in iRichards et al. I (j2006al ). Mj( z=2 ), and another where the 
value of a Q is used, Mj( Qo ), we can see in Figure 9 that there is a trend with redshift, due 
to the suggested evolution of a Q with redshift. When using a spectral index computed from 
the quasar photometry, the quasars at z ~ 2 are brighter on average with respect to the 
luminosity computed with a = —0.5, while the quasars at z ~ 4 are fainter. The difference 
in computed luminosities in the SQIID sample can be as much as ±1™8. 
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4.5. The Quasar Luminosity Function 

If there is evolution in a with redshift or luminosity, this will have a direct effect 
on the evolution of the QLF. For example, the DR3QLF is shown in Figure 10 (open 



symbo ls). If we recompute the absolute magnitudes of the bins (see 



Richards et al. 



2006al . Table 6) by using the mean redshift of the quasars in the bin in our Eq. 9 to 
compute an average a a at this mean redshift, the absolute magnitudes will shift by an 
amount AMj = 2.5(a Q + 0.5) log 10 (l + z). This term will be zero where a Q = —0.5, which 
corresponds to z = 3.1. For redshifts significantly below this, the quasar SED is steeper 
on average, and Mj is correspondingly brighter. For redshifts above 3.1, the SED's are on 
average flatter, and the Mj are fainter. 

This is shown graphically in Figure 10 where the DR3QLF is given along with points 
shifted in Mj as described above. While the space density remains essentially unchanged 
near the peak of quasar activity at z ~ 3, the points shift towards brighter Mj at lower 
redshift and towards fainter Mj at higher redshifts. For cummulative space densities as a 
function of redshift, this would mean that at lower redshifts, more objects would have Mj 
brighter than a given cutoff, increasing space densities, while at z ~ 4, fewer objects would 
make the cutoff, resulting in lower space densities. This would change the shape of the form 
in the number of quasars over time, giving rise to steeper growth at early times, with a 
more gradual decline locally. 



5. Discussion 



It has long been recognized that adopting an average spectral index for the power- 



law form of quasar spectra 



(e.g. 



Giallongo fc Vagnetti 



energy distributions can effect the evolution in the 



1992 



Francis 



1996 



Kennefick et al. 



1997 



QLF 



Wisotzki 



1998 
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Richards et al. Il2006al ). However, the adoption of an average spectral index persists in 
most studies of the QLF, with a = —0.5 being the most common value used. Recent work 
in the X-ray region has led to a general c onsensus that there is a dependence of quasar 



SED on luminosity (see 



Tang et al. 



20071 . for an exception). As for dependence on z, some 



groups find no evid ence for the evolution of X-ray spectral indices, a ox , with redshift (e.g. 



Steffen et al. 



indices with z flKelly et al. 



200 61), while others do see a linear dependence of optical to X-ray spectral 



20071 ). X-ray emission in AGN is generally taken to arise from 
a hot corona of optically thin gas heated by Compton scattering of thermal photons from 
a thick accretion disk, the likely source of the optical/UV emission. While emission from 



these two regions is likely rela ted, there is 



Kelly et al. 



i ttle e vidence for a direct link between the X-ray 



(120071 ) find no evidence for a correlation between 



and optical/UV emission, and 
auv an d a x- 

Pre v ious a ttempts to study the optical spectral energy distributions of quasars include 



Francis 



(119961 ) who fi nd a = —0.46 ± .30 fr om a sample of LBQS quasars using optical 



and NIR photometry, ICristiani fc Vio 



(119900 who compute K corrections function 



of redshift in UBV and find a 
spectra, 



Vanden Berk et al. 



-0.7 at 1000 — 5500A from their composite quasar 



(120011 ) who construct a composite quasar sp ectrum from an 



SDSS quasar sample and find a = —0.46 for the region Ly-a to H/3, and 



Pentericci et al. 



( 120031 ) who find (a) = —0.57 ± 0.33 from optical and NIR photometry of 45 z > 3.60 



SDSS quasars. There is considerable spread in the values of a within each sample, and the 



distribution of a found here and by the groups using optical and NIR photometry ( iFrancis 



1996 



Pentericci et al. 



20031 ) are remarkably similar, each having a peak around -0.5 to -0.3 



with a tail to the red that shifts the mean of the samp les to steeper values. This could 



be consistent with the findings of 



Webster et al. 



(119951 ) in their comparison of optical to 



NIR colors of radio quiet and radio loud quasars, who predict quasars should have an SED 
with a = —0.3 and that the distribution is caused by dust reddening in the host galaxies. 
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important than the evo 



ution of a a with z is the distribution of the spectral 



Perhaps more 

index values. lLa Franca fc Cristiani I (119971 ) have pointed out that, if there is a spread in 
the spectral slope, there will be a corresponding slower luminosity evolution and steeper 
QLF's. 

In gen eral, the most desirable way to present the QLF is in terms of the bolometric 



luminosity. 



Richards et al. 



(l2006bl ) demonstrate that computing bolometric luminosities 



from opt ical luminositie s assum ing a single mean quasar SED can lead to errors as large 



as 50%. iHopkins et al. 



(120071 ) have determined the bolometric QLF by combining the 
results of over two dozen quasar surveys from the hard X-ray to the mid-IR. They construct 
a model SED but allow for the distribution in the power-law components of the model, 
stressing that t here is no "effective m ean" SED. They also adopt the luminosity dependent 



value of a„ r of 



Steffen et al. 



(120061 ) . Howeve r, they do n o t ado pt a value of a ox that 



Kelly et al. 



(120071 ) and is supported by our 



depends on redshift as has been suggested by 
findings in the restframe optical. 

Given the brightness of our sample, we have not attempted to correct for the 
contribution of the host galaxy to the SED. The optical and NIR data were taken several 
years apart, and we have not considered variability in this sample. We have obtained NIR 
observations for ~ 100 more SDSS quasars in these redshift ranges, along with nearly 
simultaneous (~ 1 month) optical imaging with the aim o f adding to our survey sample and 



addressing the issue of variability. 



Richards et al. 



(120031 ) have stressed the need to correct 



for redshift dependent color effects when computing photometric spectral indexes, as failing 
to do so can lead to effects systematic with redshift. We have neglected the presence of 
emission lines in our photometric passbands. However, due to the nature of the sample - 
each quasar is sampled at essentially the same points in their SED's as is clearly seen in 
Figure 5 - the apparent evolution would persist. 
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The generally accepted view of quasar activity is ascribed to the release of gravitational 



role. The UV/optical flu x is seen 



Koratkar & Blaes 



1999, for a 



energy by accretion of matter on to a supermassive black 
as arising from a thin, optically thick accretion disk (see 
review), so a correlation between a Q and z implies evolution in the accretion process. More 
studies of the form and possible evolution of quasar SED's are obviously still needed to 
constrain theoretical models of AGN structure an d energy production m echanisms. While 



Schneider et al. 



20071 ). the form of the 



ever larger samples of quasars are discovered (e.g. 
QLF cannot be fully characterized until we have a better understanding of quasar energy 
distributions and how they are affected by luminosity, redshift, and environment. Since we 
are now coming to understand how quasar activity might be related to galaxy formation 
and evolution, quantifying the shape of the QLF and its evolution remains a pressing 
problem. Here we have presented some evidence that quasar SED's evolve with cosmic time 
and have shown that this has a direct effect on the evolution of their luminosity function. 
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Fig. 1. — SQIID JHK S filters along with the Johnson B filter "redshifted" to show where 
the restframe -B-band flux of quasars at the given redshift would be with respect to the 
corresponding SQIID NIR filter. 
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Fig. 2. — SQIID JHK S magnitudes for the sub-sample of SDSS quasars given in Table 1 
(blue triangles) along with their associated errors. Also shown are the J HK X magnitudes 



of th e 6192 quasars from the SDSS Third Data Release Quasar Catalog (ISchneider et al. 



20051 ) with measurable 2MASS detections (black dots). The SQIID data reach about l m to 



2 m fainter than the 2MASS data. 
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Fig. 3. — SDSS optical and SQIID JHK S colors for the sub-sample of SDSS quasars observed 
with SQIID given in Table 1 (blue triangles) along with their associated errors. Also shown 
are the colors of the 6192 quasars from the SDSS Third Data Release Quasar Catalog with 
measurable 2MASS detections (black dots). The dashed line shows the expected colors of 
quasars with spectral indices of a = —0.5, including emission lines and Lyman-a forest 
continuum depression. 
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magnitudes have been converted to the AB system (see text; the SDSS magnitudes are 
based on the AB system), and all magnitudes shifted to the restframe of the quasar. The 
magnitudes are also normalized to 20.0 in the bluest band completely redward of Lyman-a, 
g for the (z) = 1.88 sample, r for the (z) = 2.82 sample, and i for the (z) = 4.03 sample. 
A power law is fit to each unnormalized quasar SED separately and given in Table 2. The 
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Fig. 5. — Absolute magnitude Mj( ao ) vs. redshift for the SQIID sample (filled black symbols). 
Crosses represent quasars at z ~ 1.88, triangles at z ~ 2.82 and squares at z ~ 4.03. The 
red dots respresent SDSS quasars with 2MASS detections. 
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Fig. 6. — a) Spectral index a Q vs. redshift for the SDSS subsample observed with SQIID 
(closed symbols.) Also shown is the mean a Q for each redshift bin (open circles) with the 
value given below. The line is a linear fit to the sample given by the relation given at upper 
left. The mean for the sample as a whole is (a Q ) = —0.55 ± 0.42. b) Spectral index a Q vs. 
Mj( Qo ) for the SQIID sample. The quasars at z ~ 1.88, 2.82, and 4.03 are shown as crosses, 
triangles, and squares, respectively. The fit to the data as a function of M i(Qo ) is given in 
the upper left and shown as a solid line. 
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Fig. 7. — (Top panels) Residuals in a Q of the linear expressions given in Equations 9 and 10. 
(Bottom panels) Two dimensional projections of the residuals in a Q of the expression given 
in Equation 11. The residuals average to ~ 0, but exhibit strong scatter, with residuals as 
large as ~ 1 in all cases. 
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Fig. 8. — Left: (top) The distribution of restframe optical spectral indices, a Q , for the SDSS 
quasars detected by 2MASS. The mean of the sample is (a) = —0.70 ± 0.53 and the median 
is —0.57. (bottom) The distribution of a Q for the SQIID sample. The mean of the sample is 
(a) = —0.55 ± 0.42 and the median is —0.47. Right: The distribution of a Q for the SQIID 
sample in each redshift bin. See Table 3 for sample statistics. 
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Fig. 9. — a) The difference in absolute magnitude, Mj, assuming a = —0.5 and the a Q from 
Table 2, as a function of redshift, for the SQIID sample. Also given are the average difference 
in magnitude for each redshift interval. There is a trend with redshift reflecting the trend 
toward flatter contiuum slopes at higher redshifts for the sample, as demonstrated in Figure 
6a. b) The same difference as a function of Mj(a ). The quasars at z ~ 1.88, 2.82, and 4.03 
are shown as crosses, triangles, and squares, respectively. 
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Fig. 10.— The QLF of 



Richards et al. 



(j2006al ) (black open symbols) and the same QLF but 



with the Mi of the redshift bins calculated using an evolving a as given by Equation 9 and 
using the mean redshift of the quasars in the bin (colored filled symbols). Lower redshift 
quasars are found to be brighter and high redshift quasars are found to be fainter. This 
results in a steeper form for the evolution in quasar space densities. 
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SDSS 


094745 


.26-004113. 


.2 


2 


.8287 


18 


.922+0 


.017 


17. 


480+0. 


124 


16. 


883+0. 


202 


16. 


.270+0. 


109 


SDSS 


100423 


.27-004042. 


.9 


2 


.7320 


18. 


627+0 


.014 


17 


637+0 


067 


17. 


245+0 


.142 


16 


.869+0. 


103 


SDSS 


102832 


.09-004607. 





2 


.8592 


17 


.839+0 


.017 


16. 


986+0. 


041 


16. 


842+0. 


109 


16 


.560+0. 


.068 


SDSS 


105808 


.47+003930 


.5 


2 


.8149 


18. 


392+0 


.023 


17. 


568+0 


.065 


17. 


161+0 


.145 


16 


.342+0. 


073 


SDSS 


121323 


.94+010414. 


.7 


2 


.8292 


20. 


.185+0 


.038 


19. 


229+0. 


.203 


17. 


910+0. 


.245 


18 


.329+0. 


345 


SDSS 


121920 


.26+010736 


.1 


2 


.8005 


19 


453+0 


.022 


18. 


047+0. 


083 


17. 


323+0 


128 


17 


137+0. 


116 


SDSS 


121933 


.25+003226. 


.4 


2 


.8791 


19. 


342+0 


.030 


18. 


023+0. 


079 


17. 


394+0. 


104 


16 


611+0. 


095 


SDSS 


122730 


.37-010446. 


.1 


2 


.8701 


18. 


.798+0 


.026 


17 


780+0. 


.063 


17. 


554+0 


.114 


16 


954+0. 


105 


SDSS 


124551 


.44+010505. 


.0 


2 


.8088 


17 


896+0 


.014 


16. 


397+0 


.252 


15. 


368+0. 


153 


15 


321+0. 


184 


SDSS 


125241 


.55-002040. 


.6 


2 


.8909 


18 


.520+0 


.013 


17. 


320+0. 


150 


17. 


023+0. 


129 


16 


.492+0. 


171 


SDSS 


131128 


.35+004929 


,7 


2 


.8090 


18. 


726+0 


.015 


17. 


395+0 


059 


17. 


408+0 


120 


16 


757+0. 


062 


SDSS 


133647 


.14-004857. 


1 


2 


.7997 


17. 


.413+0 


.020 


16. 


091+0. 


049 


15. 


600+0. 


064 


15 


.360+0. 


047 


SDSS 


143307 


.40+003319 


.0 


2 


.7432 


19 


.176+0 


.022 


17. 


571+0. 


.082 


16. 


877+0. 


.131 


16 


.587+0. 


.094 


SDSS 


145754 


.03+003639. 





2 


.7603 


18. 


820+0 


.019 


17. 


810+0. 


092 


17. 


027+0. 


.134 


16 


840+0. 


.105 
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Table 1 — Continued 



SDSS DR3 designation 




Redshift 




i* 






J 






H 






K s 




SDSS 


150611 


23+001823 


6 


2 


8377 


18 


832+0 


016 


17 


466+0 


053 


17 


145+0 


119 


17 


209+0 


121 


SDSS 


094822 


96+005554 


4 


3 


8777 


20 


083+0 


037 


18 


605+0 


183 


18 


203+0 


362 


17 


559+0 


150 


SDSS 


104837 


40-002813 


6 


3 


9918 


19 


094+0 


068 


16 


219+0 


030 


15 


672+0 


049 


15 


265+0 


033 


SDSS 


105254 


59-000625 


8 


4 


1619 


19 


535+0 


024 


18 


823+0 


153 


18 


824+0 


643 


18 


255+0 


309 


SDSS 


105602 


37+003222 





4 


0361 


19 


411+0 


028 


18 


409+0 


176 


18 


160+0 


324 


17 


349+0 


140 


SDSS 


105902 


73+010404 





4 


0978 


19 


208+0 


021 


18 


122+0 


123 


17 


541+0 


200 


16 


984+0 


108 


SDSS 


110813 


85-005944 


5 


4 


0175 


19 


254+0 


020 


18 


523+0 


091 


17 


571+0 


126 


16 


717+0 


077 


SDSS 


111224 


18+004630 


3 


4 


0346 


19 


647+0 


026 


18 


667+0 


170 


18 


184+0 


361 


17 


519+0 


160 


SDSS 


120138 


56+010336 


2 


3 


8475 


19 


867+0 


028 


18 


723+0 


172 


18 


781+0 


491 


18 


291+0 


263 


SDSS 


121531 


55-004900 


4 


3 


8842 


19 


863+0 


037 


18 


660+0 


150 


17 


975+0 


219 


17 


058+0 


107 


SDSS 


122600 


68+005923 


6 


4 


2586 


18 


859+0 


021 


17 


836+0 


064 


17 


044+0 


111 


16 


490+0 


073 


SDSS 


131052 


50-005533 


2 


1 


1585 


18 


847+0 


023 


17 


650+0 


129 


17 


193+0 


225 


16 


569+0 


084 


SDSS 


135828 


74+005811 


3 


3 


9225 


19 


318+0 


022 


18 


048+0 


079 


17 


368+0 


132 


17 


240+0 


096 


SDSS 


141315 


36+000032 


3 


1 


0760 


19 


700+0 


027 


18 


352+0 


092 


17 


562+0 


186 


17 


188+0 


106 



No te, — Table[T]Optical i* photometry and redshift data are taken from the SDSS DR3Q (Schneider ct al. 
20051) . 
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Table 2. SQIID Sample Spectral Indices 



SDSS DR3 designation 




Redshift 




4 

111 




Z = 0) 




z=2) 


OLo 


&a 




.< 
«») 


SDSS 


095048 


.48-000017. 


,7 


1 


.8802 





298 


-25. 


.777 


-26. 


.068 


-1.77 





.12 


-27. 


527 


SDSS 


095938 


.28-003500. 


.8 


1 


.8753 


0. 


183 


-26 


694 


-26. 


.983 


-0.34 


0. 


09 


-26. 


798 


SDSS 


101119 


.94-004145. 


3 


1 


.8879 


0. 


234 


-26. 


225 


-26. 


.519 


-0.88 


0. 


07 


-26. 


952 


SDSS 


102517 


.58+003422. 


.0 


1 


.8879 


0. 


257 


-27 


.241 


-27 


.535 


-0.28 





05 


-27. 


281 


SDSS 


103204 


.74-001119. 


.1 


1 


.8715 





330 


-26 


.651 


-26. 


.940 


-0.39 





07 


-26. 


808 


SDSS 


103427 


.57-002233. 


.9 


1 


.8698 


0. 


375 


-26. 


205 


-26. 


.493 


-0.59 





10 


-26. 


596 


SDSS 


110725 


.70+003353. 


.8 


1 


.8732 


0. 


181 


-26 


.702 


-26 


.990 


-0.61 





.07 


-27. 


111 


SDSS 


115115 


.38+003826. 


.9 


1 


.8805 


0. 


114 


-27. 


.672 


-27. 


963 


-0.45 





.04 


-27. 


.902 


SDSS 


121655 


.39+001415. 


3 


1 


.8706 


0. 


130 


-27. 


.026 


-27 


.315 


-0.21 





.06 


-26. 


985 


SDSS 


123505 


.91—003022. 


.3 


1 


.8804 





121 


-26 


.684 


-26. 


.975 


-0.49 





05 


-26. 


967 


SDSS 


123514 


.94+004740. 


.7 


1 


.8747 


0. 


126 


-26 


263 


-26. 


.552 


-0.90 





.08 


-27. 


.007 


SDSS 


123947 


,60+002516. 


,2 


1 


.8483 


0. 


087 


-25 


.629 


-25 


.912 


-1.62 


0. 


14 


-27. 


185 


SDSS 


132742 


.92+003532. 


.6 


1 


.8736 


0. 


134 


-26 


.928 


-27. 


.216 


-0.90 





.05 


-27. 


677 


SDSS 


135605 


.41—010024. 


.4 


1 


.8860 





.264 


-26 


521 


-26. 


.815 


-0.75 





.06 


-27. 


098 


SDSS 


141015 


.36—001418. 


.9 


1 


.8758 


0. 


222 


-26. 


.630 


-26. 


.920 


-0.87 


0. 


07 


-27. 


343 


SDSS 


143641 


.24+001558. 


9 


1 


.8659 


0. 


.208 


-26 


.884 


-27. 


.171 


-0.67 


0. 


06 


-27. 


365 


SDSS 


145838 


.04+002417. 


.9 


1 


.8847 


0. 


251 


-26 


.806 


-27. 


.098 


-0.40 





.06 


-26. 


987 


SDSS 


094745 


.26-004113 


.2 


2 


.8287 


0. 


416 


-27. 


.390 


-27. 


.807 


-1.13 


0. 


09 


-28. 


725 


SDSS 


1 00423 


27—004042 


9 


2 


.7320 


0. 


281 


-27 


.553 


-27. 


974 


-0.16 





.07 


-27. 


486 


SDSS 


102832 


.09-004607. 


.0 


2 


.8592 


0. 


276 


-28. 


.440 


-28. 


.859 


0.17 


0. 


06 


-27. 


874 


SDSS 


105808 


.47+003930. 


,0 


2 


.8149 


u. 


ZIO 


97 


898 
oZo 


9fi 
-Zo. 




nil 

-U. 1 1 


0. 


.06 


97 


£89 
OoZ 


SDSS 


121323 


,94+010414. 


.7 


2 


.8292 


0. 


171 


-26 


.028 


-26 


.445 


-0.69 





.17 


-26. 


727 


SDSS 


121920 


.26+010736. 


.1 


2 


.8005 





102 


-26 


710 


-27 


.126 


-1.05 


0. 


11 


-27. 


923 


SDSS 


121933 


.25+003226. 


.4 


2 


.8791 


0. 


130 


-26 


.893 


-27 


.314 


-0.75 





.09 


-27. 


681 


SDSS 


122730 


.37-010446. 


.1 


2 


.8701 


0. 


121 


-27 


.427 


-27. 


846 


-0.12 





09 


-27. 


293 


SDSS 


124551 


.44+010505. 


.0 


2 


.8088 


0. 


103 


-28. 


.274 


-28. 


.690 


-1.12 


0. 


.07 


-29. 


591 


SDSS 


125241 


.55-002040. 


.6 


2 


.8909 


0. 


159 


-27. 


736 


-28. 


.158 


-0.46 


0. 


.08 


-28. 


093 


SDSS 


131128 


.35+004929. 


,7 


2 


.8090 


0. 


161 


-27 


467 


-27 


.883 


-0.21 


0. 


06 


-27. 


464 


SDSS 


133647 


,14-004857. 


.1 


2 


.7997 


0. 


140 


-28. 


764 


-29. 


.181 


-0.32 


0. 


04 


-28. 


918 


SDSS 


143307 


.40+003319. 


.0 


2 


.7432 


0. 


186 


-26. 


.975 


-27 


.393 


-0.66 


0. 


.08 


-27. 


628 


SDSS 


145754 


.03+003639. 


.0 


2 


.7603 


0. 


.266 


-27 


.377 


-27. 


.794 


-0.37 





.08 


-27. 


.610 
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Table 2 — Continued 



SDSS DR3 designation Redshift A u M i(z=0 ) a M i(z=2) b a a a ao M i(ao f 



SDSS 


150611 


234 


-001823 


6 


2 


8377 





309 


-27 


444 


-27 


861 


-0.36 





08 


-27 


654 


SDSS 


094822 


964 


-005554 


4 


3 


8777 





562 


-26 


977 


-27 


445 


-0.53 





22 


-27 


488 


SDSS 


104837 


40- 


-002813 


6 


3 


9918 





208 


-27 


885 


-28 


361 


-1.16 





06 


-29 


513 


SDSS 


105254 


59- 


-000625 


8 


4 


1619 





254 


-27 


552 


-28 


056 


0.45 





18 


-26 


360 


SDSS 


105602 


374-003222 





4 


0361 





212 


-27 


593 


-28 


073 


-0.16 





17 


-27 


473 


SDSS 


105902 


734 


-010404 





4 


0978 





157 


-27 


806 


-28 


298 


-0.33 





13 


-27 


999 


SDSS 


110813 


85- 


-005944 


5 


4 


0175 





243 


-27 


753 


-28 


230 


-0.48 





10 


-28 


197 


SDSS 


111224 


184 


-004630 


3 


4 


0346 





186 


-27 


346 


-27 


825 


-0.24 





16 


-27 


360 


SDSS 


120138 


564 


-010336 


2 


3 


8475 





111 


-26 


993 


-27 


454 


-0.00 





19 


-26 


600 


SDSS 


121531 


55- 


-004900 


4 


3 


8842 





096 


-27 


012 


-27 


480 


-0.85 





15 


-28 


080 


SDSS 


122600 


684-005923 


6 


4 


2586 





123 


-28 


223 


-28 


757 


-0.47 





11 


-28 


701 


SDSS 


131052 


50- 


-005533 


2 


4 


1585 





131 


-28 


188 


-28 


692 


-0.41 





11 


-28 


539 


SDSS 


135828 


744 


-005811 


3 


3 


9225 





190 


-27 


616 


-28 


095 


-0.40 





10 


-27 


913 


SDSS 


141315 


364 


-000032 


3 


4 


0760 





219 


-27 


328 


-27 


816 


-0.67 





14 


-28 


110 



Note. — Table [2] Redshift and galactic extinction data are taken from the SDSS DR3Q 



I Schn eider et al. 



2005) 



a Absolute Mi from SDSS DR3Q, column 32, X-corrected to z = 
b Absolute Mi from 



Richards et al. 



c Absolute Mi computed using values from 
using the spectral index, a , in column 6. 



( 2006aj) who if-correct to z = 2 

) who i-f-correct to z = 2, but 



Richards et al. 



(2006a 
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Table 3. Sample Spectral Index Statistics 



Quasar Sample 


N 


<*> 




median 


mean 


<T 


SDSS w/SQIID 


45 


2.81 


-27.61 


-0.47 


-0.55 


0.42 


"J" subsamplc 


17 


1.88 


-27.15 


-0.61 


-0.71 


0.43 


"H" subsample 


15 


2.82 


-27.89 


-0.37 


-0.49 


0.40 


"K" subsample 


13 


4.03 


-27.87 


-0.41 


-0.40 


0.39 



SDSS w/2MASS 6192 0.98 -25.62 -0.57 -0.70 0.53 

SDSS w/2MASS (-28 < M % < -26.5) 1336 1.44 -27.29 -0.47 -0.54 0.37 



